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Abstract. We study a class of nonlinear kinetic Fokker-Planck type equa- 
tions modeling quantum particles which obey the Bose-Einstein and Fermi- 
Dirac statistics, respectively. We establish the existence of classical solutions 
in the perturbative regime and prove exponential convergence towards the 
equilibrium. 
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1. Introduction and main results 

In recent years the rigorous mathematical study of kinetic equations has been en- 
larged to a class of models which take into account also quantum effects, cf. [53] 
for a general overview. This so-called quantum kinetic theory can be seen as an 
attempt to incorporate certain properties of an underlying quantum systems into 
the framework of classical statistical mechanics. One might hope that these "hybrid 
models" on the one hand allow for a somewhat simpler description of the particle 
dynamics, while maintaining, on the other hand, some purely quantum mechanical 
features such as generalized statistics for Bosons and Fermions. Clearly such mod- 
els can only be justified in a semi-classical regime, respectively, in situations where 
the transport properties of the particles are mainly governed by Newtonian me- 
chanics. Indeed this point of view has already been adopted in the classical paper 
by Uehling and Uhlenbeck ^24J, in which they derived their celebrated nonlinear 
Boltzmann type equation for quantum particles. 

Following their spirit most of the quantum kinetic models studied so far invoke 
nonlinear collision operators of Boltzmann type, see, e.g., [SI HH IIZ] . For these 
kind of models, a focus of interest is on the long time behavior of their solutions, 
in particular the convergence towards steady states, which generalize the classical 
Maxwellian distribution, cf. [3 [161 [HI [20] . Very often though, the simplified case 
of a spatially homogeneous gas is considered. 

In the present work, we shall also be interested in such kind of relaxation-to- 
equilibrium phenomena, in the spatially inhomogeneous case. We shall not deal 
with a Boltzmann type equation, but rather study a nonlinear Fokker-Planck type 
model (FP). More precisely, we consider 

(1.1) dtf +p-V^{f + = divp(Vp/+p/(l + K/)), 

where, for any t > 0, f — f{t,x,p) > denotes the particle distribution on phase 
space X Rp. In what follows, the spatial domain is chosen to be fix — T"^, the 
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d-dimensional torus. This setting can be seen as a convenient and mathematically 
simpler replacement for the incorporation of confining potentials V{x), needed to 
guarantee the existence of nontrivial steady states in the whole space. In we 
set K = —1 for Fermions and k = 1 for Bosons. For k = equation (jl.ip simplifies to 
the classical linear Fokker-Planck equation (or Kramer's equation) on phase space. 
For this linear model, the convergence to equilibrium has recently been studied in 
[U[Tll[T9], using several different approaches. 

The FP type model (jl.ip has been introduced in [T3] , for classical particles obeying 
an exclusion principle. A formal derivation from a generalized Boltzmann equation 
for fermions and bosons is given in [12j . and from the Uehling-Uhlenbeck equation 
in [5T]. Different physical applications can be found in [H1[I31[TS] dealing with, both, 
the spatially homogeneous as well as the inhomogeneous case (see also [10] and the 
references therein). More recently a similar but somewhat simpler FP type model 
has been proposed in [22ll23] to describe self-gravitating particles and the formation 
of Bose-Einstein condensates in a kinetic framework. The authors consider 

(1.2) dtf+p-V^f^ divp(Vp/+p/(l + /t/)), 

where, in contrast to only the diffusive part of the equation includes a non- 

linearity. As we will see both equations however share the same steady states. In 
order to deal with both models at the same time we study from now on the following 
initial value problem 

f dtf+p-VM + anf)^ divp(Vp/+p/(l-KK/)), 

I f\t=o = Mx,p), 

with a — 1, or (T — 0, corresponding to the case and (|1.2p . respectively. We 
note that the long time behavior of these models in the spatially homogeneous 
case (and for d = 1) has been rigorously investigated quite recently in [3] via an 
entropy-dissipation approach. 

In what follows, the initial phase space distribution /o G L^{T'^ x R^) is assumed 
to be normalized according to 

(1.4) // fo{x,p)dxdp^M, 

for some given mass M > 0. This normalization is conserved by the evolution. 
Moreover in the fermionic case, i.e. k = —1, we require /o(x,p) < 1, V(a;,p) £ 
T"* X R'', as usual in the physics literature [9]. In particular the latter is needed to 
define the associated quantum mechanical entropy functional, i.e. 

H[f]-^ II f^/ + /ln/-«:(l + K/)ln(l + K/))da;dp, 

which obviously requires f{t,x,p) < 1, if k = —1. It is now straightforward to 
verify that (independent of the particular choice of a) the unique steady state of 
(jl.SP is given by 

(1.5) U ^ 



(1.3) 



exp(Mi + 0) 



where the constant 9 is used to ensure that foo satisfies the mass constraint (|1.4p . 
In the bosonic case we require 9 G M+, whereas in the fermionic case we can allow 
for ^ € R, cf. [31 [7] for more details. In the latter situation the distribution (|1.5p is 
the well known Fermi-Dirac equilibrium distribution. On the other hand, for k = 1, 
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/oo is the so-called regular Bose-Einstein distribution. Finally, if k = 0, formula 
(jl.Sp simplifies to the classical Maxwellian, i.e. 

flin _ ,-|plV2 

where log(M(27r)"''/2) = ^^^^ ^j^^t in any case the equilibrium state is inde- 
pendent of X since we have chosen our spatial domain to be T**. 

In the bosonic case there is an additional difficulty, at least for c? > 3, since one 
can not associate an arbitrary large Af > to the steady state. More precisely, the 
maximum amount of mass comprised by foo is determined via 

da; dp -. M„it < oo, 

i.e. for 6 = 0. Due to mass conservation this induces a threshold on M. This 
problem, which does not appear in dimensions c? = 1 or 2, has led to the intro- 
duction of more general bosonic steady states, where an additional (^-distribution 
(appropriately normalized) is added to f oo, cf. OH]. This singular measure can 
then be interpreted as a so-called Bose-Einstein condensate (BEC). The formation 
of a (5-measure in finite or infinite time is a task of extensive research in quantum 
kinetic theory, see, e.g., OH]- For our nonlinear model though, including such 
generalized solutions on a rigorous mathematical level seems to be out of reach so 
far and we thus have to impose 9 > in the bosonic case. (Indeed, as we shall see 
below, we also require 9 > in the fcrmionic case, although for different and rather 
technical reasons.) 

Our main task here is the description of the convergence for solutions of (|1.3[) 
towards the steady state (|1.5|) . To this end we shall conceptually follow the approach 
given in jT9l where the trend to equilibrium is studied for a wide class of kinetic 
models close to equilibrium. The difference in our case being mainly that we are 
dealing with local nonlinearities which moreover are also allowed to enter in the 
transport part of the considered equation. Mathematically speaking, the approach 
is based on the so-called hypo-coercivity property of the linearized equation (26.,27j . 

We consequently linearize the solution / of p.3p around the steady state foe in the 
form 

(1-6) f ^ foo + g^/Jj^- 

where the new unknown g{t,x,p) € M can be interpreted as a perturbation of the 
equilibrium state such that 



5\/Moo dx dp = 0. 



In (II. 6p we use the additional (time-independent) scaling factor 

Moo ■ ./oo + l^foo, 

which allows for an easier description in the functional framework given below. 
Plugging (|1.6|) into (|1.3|) straightforward calculations formally yield the following 
equation for g 

(1.7) dtg + (1 + 2aKfoo)p ■ V,.g = L{g) + Q{g) . 
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Above the linearized collision operator L is given by 

L{9) = divp (Vp {g^/Jj^) + PVoo9 ^/Jj^) 

= + g l^Voo ~ \p\ I ^ + 2Kfic 

where we use the short hand notation rjoo 1 + 2k/oo- The quadratic remainder 
Q is 

(1.9) Q{g) = — ^ (divp (p^ooff^) - CMooP • Va;(5^))- 

The main result of our work is as follows. 
Theorem 1.1. Let /o be of the form 

< /o = /oo + goy/Jj^ , 

wit/i > 0. Moreover if k ^ a = 1, i.e. t/ie bosonic case with nonlinear transport, 
assume that in addition 9 > 0* , for a certain 9* > 0. 

For k (z N with k > l+d/2 there exists an > 0, such that for all fo with \\ go < 
eq, the equation (|1.3p admits a unique solution < / G C([0, cx)); _ff''(T^ x M.^)). 
Moreover 

l^-oo"\f{t)-foo) ^^<C{eo)e-^\ 
where C(eo) and t are positive constants. 

First we collect several preliminary results in the Section [2] The proof of Theorem 
II. H is then given in Section [31 

Remark 1.2. 

• For Bosons it is crucial to avoid the possible formation of a BEC by imposing 
6* > since formal calculations given in [23] indicate that an analogous 
theorem can not hold if M > Merit , see also [5] . In the fermionic case the 
reason to impose 61 > is to guarantee r]ao > 0, which we will make use of 
several times. It might be possible to overcome this restriction for fermions 
by using a different approach, see j3]. 

• The additional requirement 9 > 9* is not needed for the model ()1.2p where 
only a nonlinear diffusion operator is present. The reason for this constraint 
when dealing with (jl.ip is that we have to maintain a fundamental regu- 
larizing property in x of the transport part, of. Section [3] for more details. 
Since 9* is then determined by a transcendental equation, we do not give an 
exact value for 9* but only perform numerical experiments which indicate 
that 9* « 0.451. 

• As already discussed in [19] one could also take into account self-consistent 
potentials, which stem from a coupling to Poisson's equation. The latter 
case might be particularly interesting in semiconductor modeling, where 
the fermionic FP type equation could be used to describe the dynamical 
behavior of charge carries obeying the "physically correct" equilibrium sta- 
tistics. 

Corollary 1.3. Under the same assumptions as above 

H[f{t)] - H[f^] < Ce-^\ 
i.e. we have exponential decay in relative entropy. 
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Proof. Inserting / — /oo + ^/Moo.'? into H[f], having in mind J J ^/Jl^g dx dp = 0, we 
perform a Taylor expansion around the steady state /oo, and finally use Theorem 
11.11 This yields the assertion of the corollary . □ 



2. Study of the linearized collision operator 



We shall now derive several properties of the linearized collision operator L to be 
used in the proof of the main result. First note that L is self adjoint on L^(Rp) 
and, by partial integration, one obtains 

2 



(2-1) {L{9).9)L^{mcL) = - 



P 

^p9 + 2 "noog 



dp = — 



Moo dp . 



^''"p) onto this kernel via 



Thus the kernel of the non-positive operator L is given by 

Ker(L) = span{^/71^}. 
Let us define the orthogonal projection in L 

n(/) := (- 

where we set 

Poo = (^J Moo dp^ > 

for reasons of normalization. Note that this is only a projection in the momentum 
variable p gM.'^. Motivated by (|2.ip we introduce the following weighted space 



/Mo 



{/eL2(M^): l/L <oo}, 



where 



ll/IlL Wpff 



LI 



Wpvoo /II 



LI ' 



Here, and in what follows, we write = L'^(M.p) for simplicity. Moreover we denote 

by 



Wfh 



L^n) 



the induced norm on phase space. Obviously the Ap-norm controls the Lp-norm for 
K nonnegative. In the fermionic case (k = —1) however this is not true in general 
since rjoo may change sign. For our functional approach the control of via Ap 
is crucial and thus we have to guarantee that 7700 > by assumption. This implies 
that for K = —1 we need to impose /oo < 1/2, Vp £ M'*, or equivalently 9 > 0. This 
certainly is more restrictive than the usual bound i.e. /oo < 1, used in the physics 
literature. In summary we require > in the bosonic case to prevent BEG and 
in the fermionic case to ensure that rjoa is globally bounded away from zero. 

First we obtain a Poincare inequality for the steady state of the linearized model. 

Lemma 2.1. If 9 > then the (strictly positive and normalized) measure Moo/Poo 
satisfies a Poincare inequality on M^, i.e. 

2 



5V0 



dp 

P 



— ( 5Moodp) <Cp |Vp5pMoodp, Cp > 0. 

'00 \JR<1 / jR^i 



Note that in the fermionic case, k = — 1, this lemma indeed holds more generally 
for any 9 eR. 
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Proof. Let A be defined by ^ = — In /loo- Inequality (3.6) in [IJ shows that if A is 
uniformly convex then a Poincare inequality with measure /ioo/poo (where we divide 
by Poo for normalization purposes) holds. Moreover Theorem 3.3. in [1] allows to 
take into account also L°° perturbations of A. Thus we wish to decompose 

A = Ai +A2, 

where Ai{p) is uniformly convex and A2{p) is a L°° perturbation. To this end note 
that A is given by 

-4= -logU™— =^ + 0-21og 



(elplV2+e _ k)2 ; 2 '=UIPIV2+0_K 
' ^ ' 21og 1 



2 KelplV2+e 

We now pick Ai = |pp/2 and A2 to be the rest of the terms appearing on the right 
hand side. Then A2 is bounded if k = — 1, or if k = 1 and > which concludes 
the proof. □ 



With the above lemma in hand we can now establish the coercivity of the linearized 
collision operator. From now on we will denote by K, Ki, . . . generic (positive) con- 
stants, to be used several times in different proofs, whereas the constants Ci, C2, . . . 
will be fixed throughout the work in order to be able to track their appearance. 

Lemma 2.2. For 9 > there exists a X > such that 

{L{g),g)L2<-X\\g-Il{g)\\l^, V^GA,. 

The coercivity property (in p G M.^) of the operator L is indeed an essential require- 
ment to establish our main result. 



Proof. We start with 

{L{9),9)l2 = - Poo 

P 

< -C' 



Vr. 



/Poo 



poo J 

dp 

Poo 



2^ 1 

g dp 

Poo 



^Jj^gdp 



J ^g-}/F^ J ,/jr^gdj^ dp= -i^i 1I5 - n(5)||^2 , 



for some Ki > 0, where we used the fact that the measure Poo/Poo satisfies a 
Poincare inequality due to the previous lemma. Now to improve on the L^-norm 
we use 

(L(g),5)^. J (yp{g^n{g)) + ^rj^ig~U{g))y dp 



< -i^2||5~n(5)||^ +if3||.g-n 



Adding the two inequalities above (multiplied by appropriate constants) finishes 
the proof. □ 



We get a similar result for the derivatives w.r.t. p E M''. 
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Lemma 2.3. Let 9 > 0, then there exist positive constants Ci and C2, such that 
for any g ^ Lp with "S/ pg £ Ap 

(VpL(5), Vp5>i. < -Ci l^Ml, + C2 \\g\\ ■ 
Proof. A lengthy calculation yields 

(VpL(g), Vpff)^, = J -{Apgf + jVp^l^ (^^,7^ " ^ - 2«;Moobl') dp 
+ / TT (d^i^t^oo + 7: + 4Kd/ioo ) dp 



-^J ff^blVoo 10)7700 - 2|pp/ioo'7oo) dp. 

The last integral on the r.h.s. is dominated by the ip-norm, since /ioo decays 
exponentially fast as \p\ 00. We also have that 1/4 + 2K/ioo > C* > 0, Vp S M"*. 
This obviously holds true for the bosonic case but is also guaranteed in the fermionic 
situation where /oo < 1/2. Thus we can estimate 

{VpL{g),Vpg)^, < -Ci iVpffl^^ + -j Voo\'^ pg\^ dp + C2 \\g\l. 

and a classical interpolation argument applied to the second term on the r.h.s. 
yields the assertion of the lemma (with different constants Ci, C2). □ 

Finally we need the following technical lemma. 
Lemma 2.4. For g,h E Ap it holds that 

{m,g)L.<C,\\g\\^J\h\\^^, C3>0. 

The combination of Lemma l2.2l and Lemma |2 . 41 induces the particular choice of Ap 
and its corresponding norm. 

Proof. We first note that 

IsIIa, II^Ia, = (llVpslli. + IPVoogll^y (||Vp/i||^. + ||p?7oo/i|li.)' 
> ^ (l|Vp5lli. + /(I + bp) vlg' dp) ' [l^phWl^ + + bn vloh" dp) ' , 

for some C3 > since the Ap-norm dominates the L^-norm. Using the following 
simple algebraic estimate 

{{a^ + b^){c^ + d^))^ >ac + bd , 
(with a, 6, c, d G M) we further obtain 

C3 IsIIa, \Ma, ^ W^pdhi Wph\\L. + (^l (1 + \p\')vlg'dp 1(1 + \p\')fjlh'dp^ ' . 

The proof then follows by applying the Cauchy-Schwarz inequality to both terms 
on the right hand side and integrating by parts in the first one. □ 
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3. Convergence for the linear model and proof of Theorem 11.11 

Now we are able to establish the long time asymptotics for the linearized equa- 
tion, which eventually will be translated also to the nonlinear model (jl.Sp in the 
perturbative setting. 

Proposition 3.1. Consider the linearized Fokker Planck type equation 

(3.1) dt9+{l + 2aKU)p-V^g = L{g) , 

with L given by (II. 8p and 6 > 0. Moreover if k ~ a = 1, assume in addition that 
9 > 9*. Let the initial data go e H''{T^ x R'^), for fc e N. Then the solution g(t) 
exists globally in time and 

\\9it)'9oo\\H^ <Ce-^' , wtthC^Ci\\go\\H.),T>0, 

where the global equilibrium g^o is given by 

5oo = ( — / / go yfl^Z dp dx ) . 
We want to remind the reader of the discussion of the role of 9* in Remark 11.21 

Proof. For the proof we proceed similarly to [T^. We will sketch the main ideas 
and stress the differences which occur due to the changes in the transport operator. 
Note that, since the equation is linear, we can w.r.o.g. consider the case where 
(7oo = 0. This can always be achieved by subtracting initially the projection onto 
the global equilibrium, i.e. by considering initial data go — go — goo- 

We start with fc = 1 . The main idea of the proof is to study the time evolution of a 
combination of derivatives w.r.t. x and p. More precisely we consider the following 
functional 

T[g{t)] a \\gt + /3 || V,gf + 7 || Vp.gf + 5 (V,g, V^g) , 

where || • || denotes the standard norm on -L^(T!^ x M^) and a, /3, 7, 6 are some positive 
constants. We note that 5 has to be small enough in comparison to [3 and 7 such 
that T is positive and controlled from above and below by the square of the usual 
-ff^(T^ X M^)-norm of g. On the other hand 5 has to be strictly positive, since we 
need it in order to close the argument (see below). We aim to prove that 

(3.2) ^.^[5(i)]<-c(||5lli + l|V.,pg||y , C > 0. 

To this end we calculate the time derivatives of the various summands in T. First, 
for the L^-norm we have 

(3.3) ^\\gt ^2{L{g),g)) <-2\\\g~n{g)\l, 



where we have used that the transport part does not contribute due to its divergence 
form (which is a straightforward calculation) and the assertion of Lemma 12.21 

Next the spatial derivatives evolve according to 

(3.4) ||v,gf = 2 (V,L(g), V ^g) < -2A ||V,g - n(V,g)||^ , 



where we again used that the contribution from the transport term vanishes and 
the fact that L commutes with Vx, thus allowing us to apply Lemma 12.21 also on 
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For the derivatives w.r.t. p we get some additional terms by the coefRcient of the 
transport operator 



<2{VpL{g),Vpg)+K J J |Vpg||V,g| dpda 



Here we have used that foo as weU as Ipp/Xoo are uniformly bounded in L°°. Now 
the first term on the r.h.s. is estimated by Lemma l2.31 which yields some damping 
for Vpf; in the A-norm, i.e. 

^ l|Vp5ll' < - 2Ci ||Vp5|Ia + 2C2 \\gf + K jj |Vp5l|V,g| dpdx 
<-Ki \\Vpg\\ + K2 \V^gf + 2C2 \\gf , 
with Ki, K2 > 0. In order to deal with the term cx Igl^, we split 

g = (5-n(g))+n(g) 

and estimate 

II5II' < iff - n(g)f + ||n(5)f < \\g - n{g)f + Ct iV.gf , Ct > 0. 

In the second step, we used the classical Poincare inequality w.r. t X e T'' and the 
fact that 11(3) has zero mean on the torus, since JJ 1700 dp da; — 0. Together with 
p.Sp . this yields 

(3.6) ^ ||Vp5l|' < -i^i ||Vp.g||^ + 2C2 \\g - n(g)f + ifg ||V,g||' , 

where — K2 + 2C2CT. Note that this is an improvement to (|3.5p . since the term 
II 5 — n((7) IP can be controlled by adjusting a and having in mind (|3.3p . 

Finally we look at the mixed derivatives w.r.t. x and p, which evolve according to 

(3.7) (V,g, Vp5) = 2 (L(V,g), Vp5)-(V,g, (l + 2(jnf^ - 2aK^o,|pn) V,<?) . 

For the first term on the right hand side we invoke Lemma l2.4[ which together with 
the Cauchy-Schwarz inequality in x implies 

(L(V,g), Vpg) < C^v W.g - n(V,.g)||^ + ||Vp<?||^ , Vj; > 0. 

The second term on the right hand side of (|3.7p . which stems from the transport 
part, generates a damping for ||Va;(7|| (which the operator L can not provide since 
it only acts in p), provided that 

(3.8) l + 2Kafo.,-2Ka^lao\p\'^ >K4>0. 
Assuming for the moment that (|3.8p is true we obtain 

^ (V,g,Vpg) <2C3r/||V,5-n(V,5)llA + 2C377~M|Vp5lll-^4||V,5ll' ■ 

In summary this estimate, together with (|3.3p . p.4p . and (|3.6p . yields 
d 



- J^[git)] < - 2{\a - C21) \\g - n{g)\l - 2(/3A - C3775) ||V,5 ~ n(V,<?)||^ 

- (Xi7 - 2C^5r^) \\Vpg\l - {Ki5 - K^-f) \V^gf . 

It remains to find coefficients a,f3,j,S,ri in JF, such such that all "bad" terms in 
the above given estimates (i.e. those which come with the wrong or without sign) 
can be controlled and the differential inequality (|3.2p holds true. This can be done 
analogously to Step 4 in the Proof of Theorem 1.1 in [19] and we therefore will 
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not elaborate further on it. Note that the functional T then clearly induces a new 
norm on phase space, equivalent to i?^(R'' x T''), via || g ||^i :— J-\g\. 

To retain the (fundamental) damping property in the spatial derivatives coming 
from the evolution of the mixed term it remains to show that the constraint ()3.8|) 
holds true. We denote 

:= l + 2K/oo-2K^oob|' • 

If K = — 1 and since > 0, it is clear that ; &) > K4 > because 7700 > in 

this case. For Bosons, i.e. k = 1, however the situation is more difficult. Note that 

lim *i(p;0)= lim (?7oo -2^*00 bp) = 1, V0 > 0, 

p| — *oo \p\ — >oo 

and thus by continuity it is enough to make sure that ; 0) 7^ 0, Vp G R'', 9 > 6*. 
Straightforward calculations show that ; 6) can only be zero if 

|^elPlV2+e)'„2H2elfl'/2+^-l = , 

which implies 

Obviously this equality can not be true for 9 larger than some critical value 9*. 
Numerical experiments suggest that this critical value is approximately 9* « 0.451. 
In summary one obtains the final estimate (|3.2p . which finishes the proof for fc = 1. 

To proceed to higher order estimates in the Sobolev index A: G N, we observe that 
the proof of Lemma (|2.3p can be generalized in a straightforward way to obtain 

II l|2 2 

{dx,dp.L{g),da:,dp^g) < -Ci.k || da^^dp.glj^^ + C2.k hln';-^ > 

for any multi-indices j,l such that k ^ lil + |j| ^ 1- An induction argument in 
A: G N, similar to the one given in [TH Theorem 3.1] then yields the corresponding 
statement in . There are no additional problems due to the coefficient in the 
transport operator, since the terms containing the highest order derivatives of g 
can be treated as above and the lower order terms, which contain derivatives of 
^ K.{p ; 9), can be handled by interpolation since dp. ^'^(p ; 9) G This concludes 
the proof. □ 



Now we apply the result for the linearized equation to the nonlinear problem. 



Proof of Theorem We have to show that the quadratic nonlinearity does not 
change the estimates obtained for the linearized equation, as long as the deviation 
from the equilibrium is small. The function g — {f — /oo)/^oo solves (|1.7p from 
which we deduce 

d 2 

^ WglH" = 2 {Tg,g)-y^k + 2 {Q{g),g)j^k , 

where T := L — (1 + 2anfoo)p ■ '^x and L, Q are given in p.Sp . (|1.9p . respectively. 
From the proof of Proposition 13.11 we know that 

(Tg,g)„. < -C ( \\^-A,9 
\\j\+W<k 
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where C is, as in (13. 2|) . Thus, if we can prove the fohowing property for the nonhnear 
part 

(3.9) {Q{9),9)H-<CQ\\g\lA ^ 

\\3\ + W<k 

it follows, since || • \ — || ■ l/^t, that 



" gfnu <-'2C \ \\dx,dp^g\\'A +CQ\\gfuA ^ \dx,dp^g\\ 

\3\ + W<k I \\j\ + W<k 



At 



-^^C* Y ll^^^^p.f 111 +eC'gllffllff'«+CeC'Qll£'llL'= Yl ll^^^^Pj^il 

\\3\ + W<k ) \\j\ + W<k 

Now choosing e small enough, such that eCg < (7, we derive 

;^ II 9 1^'= < -C" I Y II^^^^K^il +C'*ll5ll?/i- Y II^^^^P.5"^ 
]\+\i\<k I \|j|+KI<fc 



This concludes the proof of Theorem II. II by maximum principle as long as Hgoll//* 
is sufficiently small. 

In order to prove p.9p . we recall that Q{g) is given by 

Q{g) = — ^ (divp(p^oo5^)) — {o^ioov ■ ^x{g^)) = Q\{g) - Qi{g), 

and note that fi^^^d^,dp^{p Hoo) e L°°{Ti x M^), for aU multi-indeces e,j € N"*. 
We shall now treat Qi{g) and ^2(9) separately, using that H''{M.'^) C L°°(M'*), 
for > d/2, together with Leibniz' formula to differentiate Q{g). It is relatively 
easy to see that p.9p holds for Qi{g), by estimates in the spirit of (j3.10|) below, 
since the A norm incorporates an additional derivate w.r.t p £ R**. The estimate 
for {Q2{g),g)uk is more complicated though, since Q2 contains a derivate w.r.t. 
X ^if^ which is not taken into account for by the A-norm. 

Thus we have to estimate 

{Q2{g),g)Hi' = -'^o- Y i^xedp. {y/JI^p ■ ^x{g^)) , d^.dp.g)^^ . 

\j\+\e\<k 

Since dx^dp. y/Hoo G L°°(T^ x M^) the highest order terms are of the form 

(VJj^P ■ ^xdx.dp^ (g^), dx,dp^9)^2 , \j\ + Kl = A: . 

Moreover, because of the additional derivative w.r.t. p in the A-norm, the most 
problematic terms are those where |^| = k. Denoting di = dx^ we compute 



{^/Jj^P ■ ^xdt{g'^), dig) ^2 = 2 {^/Jl^gp ■ Vxdeg, deg)^ 

d 



i=l 0<r<t + Si 

0<|r|<|«| + l 
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where 6i denotes the j-th standard basis vector in M''. Using divergence theorem, 
we obtain 



i^/Jj^p-^xOt^ig"^), dig)^^ = - {^/Ji^ip ■V^^g)dig, dig)^^ 



l<|r|<|« 



L2- 



To estimate the first term on the right hand side we note that (remember \t\ = k) 
(3.10) (Vm^ {p ■ Vxff) dig, dig) ^2 < Ki \\ dig ||^2 || ^xg \\ ^oo < K2\\g 1%^ \\ g l^*. , 
as soon as N 3 /c > 1 + d/2 and analogously for all other appearing terms. Since 

\j\ + W<k 

we obtain the desired estimate (13. 9p . □ 
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